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Abs tract 

An  integral  is  exhibited  which  has  the  same  value  for  all  paths  surrounding 
a  class  of  notches  in  two-dimensional  deformation  fields  of  linear  or  non-linear 
elastic  materials.  The  integral  may  be  evaluated  almost  by  inspection  for  a  few 
notch  configurations.  Also,  for  materials  of  the  elastic-plastic  type  (treated 
through  a  deformation  rather  than  incremental  formulation),  with  a  linear  re¬ 
sponse  to  small  stresses  followed  by  non-linear  yielding,  the  integral  may  be 

evaluated  in  terms  of  Irwin’s  stress  intensity  factor  when  yielding  occurs  on  a 

\ 

scale  small  in  comparison  to  notch  size.  On  the  other  hand,  the  integral  may  be 
expressed  in  terms  of  the  concentrated  deformation  field  in  the  vicinity  of  the 
notch  tip.  This  implies  that  some  information  on  strain  concentrations  is 
obtainable  without  recorse  to  detailed  non-linear  analyses.  Such  an  approach  is 
exploited  here.  Applications  are  made  to:  1)  Approximate  estimates  of  strain 
concentrations  at  smooth  ended  notch  tips  in  elastic  and  elastic-plastic  materials, 
2)  A  general  solution  for  crack  tip  separation  in  the  Barenblatt-Dugdale  crack 
model,  leading  to  a  proof  of  the  identity  of  the  Griffith  theory  and  Barenblatt 
cohesive  theory  for  elastic  brittle  fracture  and  to  the  inclusion  of  strain 
hardening  behavior  in  the  Dugdale  model  for  plane  stress  yielding,  and  3)  An 
approximate  perfectly  plastic  plane  strain  analysis,  based  on  the  slip  line 
theory,  of  contained  plastic  deformation  at  a  crack  tip  and  of  crack  blunting. 
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Introduction 

Consider  a  homogeneous  body  of  linear  or  non-linear  elastic  material 
free  of  body  forces  and  subjected  to  a  two  dimensional  deformation  field 
(plane  strain,  generalized  plane  stress,  anti-plane  strain)  so  that  all 
stresses  depend  only  on  two  cartesian  coordinates  x^  (=  x)  and 

x2  (=  y)  •  Suppose  the  body  contains  a  notch  of  the  type  shown  in  fig.  1, 
having  flat  surfaces  parallel  to  the  x  axis  and  a  rounded  tip  denoted  by  the 
arc  I*  .  A  straight  crack  is  a  limiting  case.  Define  the  strain  energy 
density  W  by 


W  =  W(x,y )  =  W(e) 


0 . > de . #  , 

13  13 


(1) 


where  e  =  [e..]  is  the  infinitesimal  strain  tensor.  Now  consider  the 
13 

integral  J  defined  by 


(2) 


Here  T  is  a  curve  surrounding  the  notch  tip,  the  integral  being  evaluated 
in  a  contraclockwise  sense  starting  from  the  lower  flat  notch  surface  and 
continuing  along  the  path  T  to  the  upper  flat  surface.  X  is  the  traction 
vector  defined  according  to  the  outward  normal  along  T  ,  =  °^jnj  >  u 

is  the  displacement  vector,  and  ds  is  an  element  of  arc  length  along  T  . 

We  now  prove  that  the  integral  J  has  the  same  value  for  all  paths  T  . 
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This  path  independence,  combined  with  the  fact  that  J  can  often  be  directly 

evaluated,  is  the  key  to  a  variety  of  subsequent  applications  to  non-linear 

strain  concentration  problems.  To  prove  path  independence,  consider  any 
&  & 

closed  curve  T  enclosing  an  area  A  in  a  two-dimensional  deformation 
field  free  of  body  forces.  An  application  of  Green's  Theorem  [1]  in  cartesian 
coordinates  x  =  x  ,  x9  =  y  gives 


3W  _  3 _ 

3x  3x. 
3 


3u. 

sr03  dxdy  • 


(3a) 


Differentiating  the  strain  energy  density. 
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3e.  .  3e . . 

=  o..  — il 
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ij  3x 


(by  eq.  1) 


1  a  a  3u-s 

k  a  [L_  ( — i.)  +  (  .  .1 )  ] 

2  ijL3x  v3x.  3x  ,‘3x. 


.  3u. 

_  3 _  ( _ lx 

ij  3x .  v3x  1 
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(since  o. .  =  o . . ) 
13  31 


a  , 

dx.  loij  ax  ' 

3  J 


3a . . 

(since  ~  °) 
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(3b) 


The  area  integral  in  eq. 


*  (Wdy  -  T  • 


3u, 

3x 


(3a)  vanishes  identically,  and  thus 


ds)  =  0  for  any  closed  curve  f 


(3c) 


Path  independence  of  the  integral  J  defined  by  eq.  (2)  is  now  a  straight¬ 
forward  consequence.  Consider  any  two  paths  and  surrounding  the 
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notch  tip,  as  does  T  in  fig.  1.  Traverse  in  the  contraclockwise 

sense,  continue  along  the  upper  flat  notch  surface  to  where  intersects 

the  notch,  traverse  in  the  clockwise  sense,  and  then  continue  along  the 

lower  flat  notch  surface  to  the  starting  point  where  intersects  the  notch. 

This  describes  a  closed  contour  so  that  the  integral  of  Wdy  -  T  •  (3u/3x)  ds 
vanishes.  But  T  =  0  and  dy  =  0  on  the  portions  of  path  along  the  flat  notch 
surfaces.  Thus  the  integral  along  contraclockwise  and  the  integral  along 

clockwise  sum  to  zero.  J  has  the  same  value  when  computed  by  integrating 
along  either  or  T ,  and  path  independence  is  proven.  We  assume,  of 

course,  that  the  area  between  curves  and  is  free  of  singularities. 

Our  scheme  of  approach  to  notch  strain  concentrations  is  first  to  point 
out  in  the  next  section  that  J  can  be  explicitly  evaluated  for  a  variety  of 
configurations  and  cases.  The  relation  of  the  integral  J  to  the  rate  of 
change  of  potential  energy  with  respect  to  notch  length  will  also  be  noted. 

Then,  in  subsequent  sections,  by  appropriate  choices  of  the  curve  T  on  which 
J  is  evaluated,  a  variety  of  notch  deformation  concentration  problems  will  be 
discussed. 

Evaluations  of  the  Integral  J 

Two  special  configurations.  The  integral  J  may  be  evaluated  almost 
by  inspection  for  the  two  special  configurations  shown  in  fig.  2.  These 
configurations  are  not  of  great  practical  interest,  but  are  useful  in  illust¬ 
rating  the  relation  between  the  integral  J  and  potential  energy  variation 
rates  which  is  taken  up  below.  In  fig.  2a,  a  semi-infinite  flat  surfaced 
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notch  in  an  infinite  strip  of  height  h  ,  loads  are  applied  by  clamping  the 
upper  and  lower  surfaces  of  the  strip  so  that  the  displacement  vector  y 
is  constant  on  each  clamped  boundary.  Take  T  to  be  the  dashed  curve  shown 
which  stretches  out  to  x  =  *  »  .  There  is  no  contribution  to  J  from  the 
portion  of  T  along  the  clamped  boundaries  since  dy  =  0  and  3u/3x  =  0 
there.  Also  at  x  =  -  »  ,  W  =  0  and  3u/3x  =  0  .  The  entire  contribution 
to  J  comes  from  the  portion  of  T  at  x  =  +  «®  ,  and  since  3u/3x  =  0  there, 

J  =  W  h  (4) 

00 


where  is  the  constant  strain  energy  density  at  x  =  +  »  . 

Now  consider  the  similar  configuration  of  fig.  2b,  with  loads  applied 
by  couples  M  per  unit  thickness  on  the  beam  like  arms  so  a  state  of  pure 
bending  (all  in-plane  stresses  vanishing  except  oxx)  results  at  x  =  -  »  . 
For  the  contour  T  shown  by  the  dashed  line,  no  contribution  to  J  occurs 
at  x  =  +  ®  as  W  and  T  vanish  there,  and  no  contribution  occurs  for  por¬ 
tions  of  T  along  the  upper  and  lower  surfaces  of  the  strip  as  dy  and  T 
vanish.  Thus  J  is  given  by  the  integral  across  the  beam  arms  at  x  =  -  00 
and  on  this  portion  of  T  ,  dy  =  -  ds  ,  T^  =  0  ,  and  Tx  =  "  °xx  •  We  end 
up  integrating 


3ux 


a  t — -  -  W  =  a  e  -W  =  a . .  e .  . 

xx  3x  xx  xx  i] 


W 


e. .da. .  =  fi 


(5a) 
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across  the  two  beam  arms,  where  ft  is  the  complimentary  energy  density. 
Thus,  letting  Q^(M)  be  the  complimentary  energy  per  unit  length  of  beam 
arm  per  unit  thickness  for  a  state  of  pure  bending  under  moment  per  unit 
thickness  M  , 


J  =  20b(M)  .  (5b) 

Our  motives  will  be  clear  enough  if  we  now  employ  path  independence  to 
shrink  the  curve  T  down  to  the  curved  notch  tip,  denoted  by  I*  in  fig.  1. 
Since  the  traction  vector  T  vanishes  on  the  notch  surface 


Without  solving  generally  non-linear  and  unmanageably  difficult  boundary  value 
problems,  we  have  extracted  a  very  significant  bit  of  information  on  the  con¬ 
centrated  deformation  field  at  the  notch  tip.  Recognizing  W  as  a  function 
of  the  surface  strain  along  the  tip,  we  obtain  an  averaged  value  of  the  con¬ 
centrated  strain  by  simply  being  able  to  estimate  the  strain  energy  in  a  homo¬ 
geneous  deformation  field,  eq.  (4),  or  to  solve  the  problem  of  pure  bending  of 
a  beam,  eq.  (5b).  Such  simple  connections  are  the  essence  of  the  method  pre¬ 
sented,  but  we  defer  the  details  to  subsequent  sections  and  take  up  next  the 
estimate  of  J  for  small  zones  of  notch  tip  yielding  in  materials  of  the 
elastic-plastic  type. 


independence  applies  only  to  elastic  materials  for  which  a  strain  energy 
density  exists.  Thus  we  are  here  forced  to  consider  deformation  theories  of 
plasticity  (non-linear  elasticity)  rather  than  the  physically  more  appropriate 
incremental  theories.  This  should  not  be  a  great  drawback  [2,3]  for  problems 
of  monotonic  loading.  Specifically,  we  shall  here  understand  the  phrase 
"elastic-plastic  material"  to  mean  an  elastic  material  exhibiting  a  linear 
Hookian  response  to  stresses  inside  some  initial  yield  surface  in  stress  space, 
and  a  non-linear  (strain  hardening)  response  to  stresses  outside  the  initial 
yield  surface.  We  consider  a  narrow  crack-like  notch  and  derive  a  formula  for 
the  integral  J  when  the  yielded  zone  is  confined  to  a  small  region  near  the 
notch  tip  that  is  negligible  in  size  in  comparison  to  geometric  dimensions  such 
as  notch  length,  unnotched  specimen  width,  etc. 

The  situation  envisioned  has  been  termed  "small  scale  yielding"  [2,4,5], 
and  we  employ  a  special  boundary  layer  type  formulation  of  the  problem  dev¬ 
eloped  in  work  by  the  author  [2,5]  and  others  [6,7],  and  described  most  fully 
in  a  recent  survey  [8]  on  crack  plasticity.  The  essential  ideas  are  illustrated 
in  reference  to  fig.  3.  Loadings  symmetrical  about  the  narrow  notch  are 
imagined  to  induce  a  deformation  state  of  plane  strain.  First,  consider  the 
linear  elastic  solution  of  the  problem  when  the  notch  is  presumed  to  be  a 
sharp  crack.  Employing  polar  coordinates  r,9  with  origin  at  the  crack  tip, 
the  form  of  stresses  in  the  vicinity  of  the  tip  are  known  [9-11]  to  exhibit 
a  characteristic  inverse  square  root  dependence  on  r  : 

KI 

a..  =  - =- jtt  f.  .(e)  +  other  terms  which  are  bounded 

«  (2i,r)1/2  ” 


at  the  crack  tip. 


(7a) 
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Here  the  set  of  functions  f„(0)  are  the  same  for  all  symmetrically  loaded 
crack  problems,  and  for  the  isotropic  material  [9,10,12] 


f  (6)  =  cos  (6/2)  [1  -  sin  (0/2)  sin  (30/2)] 

XX 

f  (0)  =  cos  (0/2)  [1  +  sin  (8/2)  sin  (30/2)]  (7b) 

yy 

fxy(0)  =  fyx(0)  =  sin  (0/2)  cos  (0/2)  cos  (30/2) 


Kj  is  Irwin's  stress  intensity  factor  [10,12],  a  function  of  geometry  and 
applied  load,  depending  linearly  on  the  latter,  that  has  been  tabulated  [12]  for 
a  wide  variety  of  crack  configurations  on  the  basis  of  extensive  recent  elastic 
crack  solutions  generated  in  studies  on  fracture  mechanics.  Now  suppose  the 
material  is  elastic-plastic,  the  flat  surfaced  notch  is  either  a  sharp  crack 
or  a  narrow  void,  and  the  load  level  is  sufficiently  small  so  that  a  yield  zone 
forms  near  the  tip  which  is  small  compared  to  notch  length  and  similar  character¬ 
istic  dimensions  (small  scale  yielding,  fig.  3a).  One  anticipates  that  the 
elastic  singularity  governs  stresses  at  distances  from  the  notch  root  that  are 
large  compared  to  yield  zone  and  root  radius  dimensions  but  still  small  com¬ 
pared  to  characteristic  geometric  dimensions  such  as  notch  length.  The  actual 
configuration  of  fig.  3a  is  then  replaced  by  the  simpler  semi-infinite  notch 
in  an  infinite  body  (fig.  3b),  and  a  boundary  layer  approach  is  employed  re¬ 
placing  actual  boundary  conditions  in  fig.  3a  with  the  asymptotic  boundary  con¬ 
ditions 


a. . 


(2irr) 


1/2 


as 


v  °° 


5 


(7c) 
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where  is  the  stress  intensity  factor  from  the  linear  elastic  crack 

solution. 


Such  boundary  layer  solutions  for  cracks  are  mathematically  exact  in 
the  plastic  region  only  to  the  first  nonvanishing  term  of  a  Taylor  expansion 
of  complete  solutions  in  the  applied  load  [2,4,5].  But  comparison  [8]  with 
available  complete  solutions  indicates  the  boundary  layer  approach  to  be  a 
highly  accurate  approximation  up  to  substantial  fractions  (typically,  one- 
half)  of  net  section  yielding  load  levels.  We  now  evaluate  the  integral  J 
(eq.  2)  from  the  boundary  layer  solution,  taking  r  to  be  a  large  circle  of 
radius  r  in  fig.  3b  : 


J 


r+ir 

r  [W(r,6>  cos  6  -  T  (r,0) 
* -v 


(r,e)]  de  . 


(8a) 


By  path  independence  we  may  let  r  -*•  «*  and  since  W  is  quadratic  in  strain 
in  the  elastic  region,  only  the  asymptotically  approached  inverse  square  root 
elastic  stress  field  (eq.  7c)  contributes.  Working  out  the  associated  plane 
strain  deformation  field,  one  finds 


J 


1  -  v2  2 
E  KI 


for  small  scale  yielding  , 


(8b) 


where  E  is  Young's  modulus  and  v  Poisson's  ratio.  We  have  thus  evaluated 
J  for  an  infinite  variety  of  problems ,  and  made  an  unexpected  use  of  linear 
elastic  crack  solutions  in  elastic-plastic  problems. 
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Primarily,  we  will  later  deal  with  one  configuration,  the  notch  of 
length  2a  in  a  remotely  uniform  stress  field  (fig.  4).  Here  [10] 

2 

„  /  .1/2  j  _  tt(  1  -  v  )  2  c  ..  . 

KT  =  a  (ira)  and  J  =  - = -  a  a  for  small  scale 

I  00  E  00 

yielding  (9) 

2 

For  plane  stress,  the  same  result  holds  for  J  with  1  -  v  replaced  by 
unity.  Thus  far  only  loadings  symmetrical  about  the  narrow  notch  line  have 
been  considered  (the  opening  mode  in  Irwin's  terminology  [10]).  The  same 
computation  may  be  carried  out  for  more  general  loadings.  Letting  , 

Kjj  ,  and  Kjjj  be  elastic  stress  intensity  factors  [10]  for  the  opening, 
in-plane  sliding,  and  anti-plane  sliding  modes,  respectively,  of  notch  tip 
deformation  one  readily  obtains 

2 

,  1  -  v  fv 2  ,  Jl  .  ,  1  +  v  „2  ,  ,,  - 

J  =  — g -  (Kj  +  Kjj)  +  — g —  (small  scale 

yielding)  (10) 


A  physical  interpretation  in  terms  of  energy  comparisons  for  notches 
of  neighboring  size.  In  view  of  the  path  independence  it  should  come  as  little 
surprise  that  a  simple  physical  interpretation  can  be  attached  to  J  .  Con¬ 
sider  the  two  dimensional  deformation  of  a  body  with  cross  section  A'  and 
bounding  curve  T'  .  Letting  T"  be  that  portion  of  I”  on  which  non¬ 
vanishing  tractions  T  are  prescribed,  we  define  the  potential  energy  per 
unit  thickness  by 


■*r 


f! 


P 


f 


W  dxdy 


T  •  u  ds 


(11) 
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Suppose  the  body  contains  a  flat  surfaced  notch,  as  in  fig.  1,  with  a  tip 
at  x  =  A  and  denote  the  potential  energy  by  P(  JZ.)  .  We  compare  this  with 
the  potential  energy  P(S.  +  LI)  of  a  body  identical  in  every  respect  (same 
shape  and  same  distribution  of  prescribed  displacements  and/or  tractions  on 
the  boundary),  except  that  the  notch  is  longer  with  the  tip  now  at  a  neigh¬ 
boring  position  x  =  £  +  A£  .  The  shape  of  the  curved  notch  tip  r  is 
imagined  the  same  in  both  cases.  Then  one  may  show  that  the  integral  J 
defined  by  eq.  (2)  is 

T  lim  pU  +  AQ  -  P(&)  3P  ,  . 

"  n+o  a z  "  '  dz  KXZ} 

the  rate  of  decrease  of  potential  energy  with  respect  to  notch  length  under 
fixed  prescribed  boundary  values.  The  proof  is  fairly  long,  and  is  deferred 
to  a  subsequent  report  both  in  the  interests  of  brevity  and  because  a  direct 
focus  on  energy  variations  is  a  subject  in  itself  allowing  a  variety  of  fur¬ 
ther  results  on  non-linear  strain  concentrations  at  notches .  For  the  present , 
we  simply  note  below  that  J  does  indeed  equal  the  potential  energy  decrease 
rate  for  the  cases  evaluated  above. 

When  loading  is  entirely  by  prescribing  boundary  displacements  as  in 
fig.  2a,  the  potential  energy  equals  the  strain  energy  and  we  check  from 
eq.  (4)  that  J  is  indeed  equal  to  the  rate  of  decrease  of  strain  energy 
with  respect  to  notch  length.  When  loading  is  entirely  by  prescribing  boundary 
tractions  as  in  fig.  2b,  a  simple  application  of  virtual  work  in  eq.  (11) 
shows  that  the  potential  energy  equals  minus  the  complimentary  energy,  and 
eq.  (5b)  checks  J  to  be  the  rate  of  increase  of  complimentary  energy  with 
respect  to  notch  length.  Readers  familiar  with  elastic  fracture  mechanics 
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will  recognize  eqs.  (8b-10)  as  giving  results  for  J  which  are  identical 
to  the  potential  energy  decrease  rate  calculated  in  general  form  by  Irwin 
[10]  in  terms  of  his  stress  intensity  factors  for  linear  elastic  crack  pro¬ 
blems.  When  deviations  from  linearity  occur  over  a  region  negligible  in  com¬ 
parison  to  characteristic  dimensions  of  a  problem  (small  scale  yielding),  it 
is  clear  that  the  energy  decrease  rate  should  be  insensitive  to  the  yielding 
and  thus  still  given  by  the  linear  elastic  expressions  so  that  we  again  check 
eq.  (12).  The  physical  interpretation  of  J  in  terms  of  a  variation  rate  in 
potential  energy  is  useful  in  enabling  an  approximate  treatment  of  notch  strain 
concentrations  when  yielding  is  not  confined  to  a  small  region  near  the  tip. 

For  example,  suppose  we  can  solve  a  simplified  elastic-plastic  model  for 
yielding  near  a  notch.  Even  though  the  solution  may  be  very  wrong  in  detail, 
we  might  confidently  expect  the  model  to  accurately  predict  a  gross  feature 
of  the  solution  such  as  the  rate  of  variation  of  potential  energy  with  respect 
to  notch  length.  Thus  we  compute  J  from  the  simplified  model,  and  then 
ignore  the  model  and  go  on  to  employ  techniques  of  the  next  sections  to 
analyze  strain  concentrations  at  the  notch  tip.  Examples  are  discussed 
subsequently . 

Strain  Concentrations  at  Smooth  Ended  Notches 

When  we  shrink  the  contour  T  down  to  the  curved  tip  of  a  smooth 

ended  notch,  as  in  deriving  eq.  (6),  there  results 

r+rr/2 

Wdy  =  W((|>)  r  (<j>)  cos  4>  d<{>  . 

rt  J-ir/2 


J 


(13) 


-13- 


Here  <|>  is  the  angle  between  the  tangent  line  at  a  point  on  the  curved 
tip  and  the  y  axis  (fig.  5)  and  r^(^>)  is  the  radius  of  curvature  at 
that  point.  Viewing  J  as  a  known  function  of  the  applied  load,  we  obtain 
an  approximate  estimate  of  the  maximum  concentrated  strain  by  choosing  some 
reasonable  dependence  of  the  strain  energy  on  <J>  .  Some  guidance  on  what 
constitutes  a  reasonable  choice  is  given  by  a  perusal  of  known  solutions 
for  notches . 

A  well  known  [13]  feature  of  the  linear  elastic  problem  of  an  ellip¬ 
soidal  inclusion  imbedded  in  a  matrix  of  another  material,  subjected  to 
a  uniform  remote  stress  state,  is  that  the  inclusion  undergoes  a  homogeneous 
deformation.  Strains  in  the  inclusion  are  independent  of  position.  The  ellip¬ 
soidal  void  is  a  limiting  case  with  an  inclusion  of  zero  elastic  modulus,  sur¬ 
face  strains  being  compatible  with  a  homogeneous  deformation  of  the  imagined 
zero-modulus  inclusion.  Now  consider  the  special  case  of  a  narrow  ellip¬ 
tical  void  symmetrically  loaded  under  plane  strain  conditions.  Let  emax 
be  the  maximum  extensional  strain,  occurring  at  the  notch  tip  on  the  semi¬ 
major  axis.  Then  the  imagined  inclusion  undergoes  a  homogeneous  strain 

e  =  e  and  e  =  e  .  ,  where  e  .  is  the  surface  extensional  strain 

yy  max  xx  min  min 

on  the  semi -minor  axis.  The  latter  is  negligible  in  comparison  to  emax  for 
a  narrow  ellipse  and  essentially  the  inclusion  undergoes  a  single  homogeneous 
strain  e  in  the  y  direction.  Matching  surface  extensional  strain  e 
with  extensional  strain  in  the  imagined  zero  modulus  inclusion  in  a  direction 
tangent  to  the  void  boundary, 

e  =  e  cos  d>  when  e  .  is  negligible  ,  (14a) 

max  T  min  s  6  * 
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where  <j>  is  the  tangent  angle  introduced  in  fig.  5. 

As  an  approximation  for  our  flat  surfaced  notch  problems,  we  may  con¬ 
sider  the  deformation  of  the  notch  surface  near  the  tip  to  be  describable 
by  the  homogeneous  straining  of  an  imagined  zero  modulus  inclusion.  Thus 
the  maximum  concentrated  strain  is  estimated  from  eqs.  (13,14a)  by  setting 

r+ir/2  2 

J  ss  W(e  cos  <f) )  r.  (4>)  cos  $  d<|>  .  (14b) 

J-tt/2  max  T 

Here  W(e)  is  the  strain  energy  accompanying  an  extensional  strain  e  on 
the  notch  surface.  One  might  instead  presume  that  the  distribution  of  strain 
energy  with  4>  in  the  linear  elastic  ellipse  problem, 

4 

W  =  W(e  )  cos  <f>  when  e  .  is  negligible  ,  (15a) 

max  T  mm  6  6  * 

is  more  likely  to  be  a  better  approximation  for  our  flat  surfaced  notch  pro¬ 
blems.  In  this  case  the  maximum  strain  is  estimated  from 

r+ir/2 

J  »  W(e  )  cos  r  (4)  d<f>  .  (15b) 

max  J-it/2  1 

Certainly,  neither  assumption  is  exact  and  there  is  no  way  (without 
recourse  to  the  detailed  solutions  we  are  trying  to  avoid)  to  decide  which 
gives  the  best  answer.  As  we  shall  see,  resulting  answers  for  maximum  strain 
are  somewhat  sensitive  to  our  choice.  An  examination  of  some  results  from 
known  anti-plane  strain  solutions  does,  however,  suggest  that  eq.  (14b)  be 
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employed.  If  we  consider  the  elliptical  void  subjected  to  symmetrical  anti¬ 
plane  strain  loadings ,  the  imagined  zero  modulus  inclusion  undergoes  the  homo¬ 
geneous  strain  e  =  y  /2  where  y  is  the  maximum  concentrated  "engin- 
yz  max  max 

eering"  shear  strain  occurring  at  the  notch  tip.  Matching  shear  strains  y 
on  the  notch  surface  with  shears  in  the  inclusion 

y  =  y  cos  (J>  .  (16a) 

max  Y 

It  happens  that  a  solution  for  a  crack  in  a  uniform  remotely  applied 
anti-plane  stress  field  has  been  given  [2,6]  for  elastic-plastic  materials 
with  arbitrary  non-decreasing  relations  between  stress  and  strain  in  the  plastic 
range.  Neuber  [6]  first  pointed  out  that  such  crack  solutions  also  give  sol¬ 
utions  for  a  family  of  smooth  ended  notches  since  stress  trajectories  in  the 
crack  solution  may  be  thought  of  as  traction  free  notch  surfaces.  Methods  of 
locating  such  notch  boundaries  are  discussed  in  [8],  If  we  take  the  small 
scale  yielding  solution  [2]  in  which  the  crack  is  viewed  as  semi-infinite  and 
generate  the  family  of  smooth  ended  semi-infinite  notches  by  methods  of  [8], 
the  surprising  result  turns  out  that  eq.  (16a)  is  satisfied  on  the  notch  sur¬ 
face  regardless  of  the  relation  between  stress  and  strain  in  the  plastic  range. 
Recalling  that  small  scale  yielding  solutions  give  details  correctly  only  near 
the  notch  tip,  the  assumption  of  a  homogeneously  deformed  zero  modulus  inclusion 
accurately  describes  deformations  of  the  notch  tip  at  least  at  low  load  levels. 
The  shape  of  notch  surfaces  generated  by  stress  trajectories  of  the  crack  sol¬ 
ution  do  depend  on  the  stress-strain  relation.  Parabolas  result  for  linear 
elastic  behavior.  Perfectly  plastic  behavior  leads  to  a  semi-circular  notch 
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tip  [8]  where  the  stress  trajectories  pass  through  the  crack  plastic  zone, 
and  relatively  flat  notch  surfaces  where  the  trajectories  pass  through 
elastic  material.  This  is  especially  encouraging  as  our  formulation  is  for 
flat  surfaced  notches  and  our  main  interest  is  with  ductile  materials. 

Neuber  [6]  pointed  out  for  this  class  of  anti-plane  strain  problems 
that  the  product  of  stress  and  strain  concentration  factors  is  independent 
of  the  stress-strain  relation.  A  more  careful  examination  of  mathematical 
results  [8]  indicated  the  Neuber  result  to  be  valid  only  for  small  scale 
yielding  of  elastic-plastic  materials,  and  only  for  the  comparison  of  stress 
and  strain  concentrations  at  notch  tips  whose  geometry  does  depend  on  the 
stress-strain  relations.  Thus  the  Neuber  result  is  only  approximately 
correct.  Rather,  our  present  considerations  of  path  independent  integrals 
and  their  relation  to  potential  energy  variations  appear  to  be  the  embracing 
principles  which  lead  to  Neuber' s  result  in  special  cases.  If  we  approximate 
notch  tip  deformations  by  eq.  (16a)  in  anti-plane  strain  problems,  the  maxi¬ 
mum  concentrated  strain  is  estimated  from 

f+it/2 

J  ^  W(y  cos  <|>)  r.  ( <t> )  cos  4>  d<f>  ,  (16b) 

J-,/2  malt  t 

where  W(y)  is  the  strain  energy  accompanying  an  anti-plane  strain  y  of 
the  notch  surface.  Strict  equality  holds  for  notch  surfaces  generated  by 
stress  trajectories  of  the  small  scale  yielding  crack  solution,  with  J  set 
equal  to  its  linear  elastic  value  (last  term  of  eq.  10). 
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Having  received  some  encouragement  from  the  anti-plane  strain  case, 
we  now  proceed  to  use  eq.  (14b)  (homogeneous  deformation  of  void  interior 
in  the  neighborhood  of  the  notch  tip)  to  estimate  concentrated  deformations 
in  the  plane  straining  of  a  flat  surfaced  notch  with  a  semi-circular  tip 
( r ^ ( <(> )  =  r  =  a  constant).  We  shall  always  first  give  results  in  terms  of 
J  and  then,  to  put  equations  in  more  familiar  terms,  give  results  for  the 
special  value  of  J  in  eq.  (9)  for  small  scale  yielding  near  the  narrow 
notch  of  length  2a  in  a  uniform  remote  stress  field  ow  (fig.  4).  The 
reader  should  recall,  however,  that  J  is  known  either  approximately  or 
exactly  for  a  large  number  of  other  cases  and  configurations. 

For  linear  elastic  behavior  stress-strain  relations  and  strain  energy 
density  on  the  notch  surface  under  plane  strain  conditions  are 

a  =  E  -  ■,  e  ,  W(e)  =  - — . ■=—  •  (17a) 

1  -  v  2(1  -  v  ) 


Inserting  in  eq.  (14b)  and  solving  for  the  maximum  concentrated  strain, 


max 


15(1  -  v  )  J 
8E  r_ 


1/2 


(1  -  v  )  o_ 


=  2.43 


(V 


1/2 


(17b) 


The  numerical  factor  of  2.43  appears  appropriate  for  a  flat  surfaced  notch 
when  one  recalls  the  factor  2.00  for  the  narrow  elliptical  hole  [14]  of  the 
same  root  radius.  The  ellipse  represents  a  much  less  abrupt  transmission  of 
load  around  the  notch  tip  (the  ratio  of  semi-minor  axis  length  to  root  radius 
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approaches  infinity  as  the  ellipse  is  narrowed  toward  zero  thickness). 

For  perfectly  plastic  behavior,  let  Oy  be  the  yield  stress  at  points 
along  the  notch  surface  under  plane  strain  conditions.  With  elastic  com¬ 
pressibility  a  limiting  yield  stress  in  plane  strain  is  really  approached 
asymptotically  when  plastic  strains  become  large  compared  to  elastic  strains 
[15],  but  we  shall  consider  a  simple  picture  of  linear  behavior  up  to  Oy  and 
constant  stress  thereafter.  Then 

o  =  Oy  ,  W(e)  =  Oy  (e  -  Ey)  +  ^  OyEy  ,  for  e  >  Ey  ,  (18a) 

and  eq.  (17a)  applies  for  e  <  Ey  ,  where 

(1  -  V2)  Oy 

ey  =  - g -  (18b) 


is  the  yield  strain  given  by  Hooke's  law  under  plane  strain  conditions.  One 
then  obtains  from  eq.  (14b),  when  the  yield  strain  is  exceeded. 


2  1/2  5/2 

G  G  G 

(_maX)  (_Y_)  ^  max 

E  £„ 

Y  max  Y 


15  J 
8°YeYrt 


15  TT 


(~) 


a_ 

r. 


(18c) 


This  equation,  along  with  the  linear  elastic  result,  is  plotted  as  the  solid 
curve  in  fig.  6.  The  ratio  of  maximum  concentrated  strain  to  yield  strain  is 
shown  as  a  function  of  the  square  root  of  the  terms  on  the  right  in  eq.  (18c), 
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since  this  dimensionless  representation  of  the  applied  load  has  the  value 

unity  at  initial  yield  and,  at  least  at  low  load  levels  when  J  is  given  by 

the  linear  elastic  expression  of  eqs.  (8,9),  is  linear  in  the  applied  load. 

The  left  side  of  the  above  equation  may  be  developed  in  a  series,  and  if  we 

neglect  all  terms  which  vanish  when  e  /ew  >>  1  , 

6  max  Y 


max 


3  ,  J 

4  ^eY  +  o,,r 


Yt 


.075 


+  2.36  e„  ( —  ) 

Y  0  y 


(18d) 


This  equation  gives  a  result  for  the  maximum  concentrated  strain  which  is  15% 
too  large  at  the  initial  yield  load;  the  discrepancy  with  eq.  (18c)  becomes 
imperceptible  at  loads  greater  than  three  times  the  initial  yield  load. 

Strain  hardening  poses  no  special  difficulty  with  our  present  method. 

For  simplicity  of  resulting  formulae,  we  shall  limit  attention  to  cases  where 
concentrated  strains  are  very  large  compared  to  the  initial  yield  strain.  Thus, 
for  linear  work  hardening  in  the  relation  between  stress  and  strain  on  the  notch 
surface , 


o  =  ay  + 


w 


1  -  v 


(e  -  ev) 


e  >  e. 


(19a) 


eq.  (14b)  leads  to 


max 


5  E  , 

4  E  EY  1 
w 


1  —  ( 
4  E  eY  1 
w 


1  + 


6  E  J 
w 


5  £  Eo^r 


1/2 


-  1) 


1  + 


c  E 

6tt  w 

5  E 


o 

(—) 


il/2 
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(19b) 
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Here  E  is  the  tangent  modulus  (effectively)  and  terms  of  order  avev 
and  Eweye  have  been  neglected  in  the  expression  for  the  energy  density,  as 
appropriate  for  large  strains  and  for  Ew  <<  E  .  For  a  power  law  relating 
stress  and  strain  on  the  notch  surface  in  the  hardening  range, 


o 


aY  (e/eY) 


N 


e  >  e. 


(20a) 


eq.  (14b)  results  in 


'max  ~  Y 


(N+l/2)  (N+3/2)  T  (N+l/2) 
r  (1/2)  T  (N+l) 


eY°Yrt 


1/ ( 1+N ) 


tt  (N+l/2)  (N+3/2)  r  (N+l/2)  (V  a__ 

°Y  rt 


r  (i/2)  r  (N+i) 


1/(1+N) 


(20b) 


where  T( ..)  is  the  Gamma  Function.  Again  terms  of  order  have  been 

neglected  in  the  energy  density. 

The  reader  is  cautioned  again  that  the  latter  forms  in  expressions  above 
for  emax  apply  only  to  the  vase  of  small  scale  yielding  near  the  narrow 
notch  of  length  2a  in  fig.  4  (that  is,  for  aa  less  than  approximately  half 
of  ©Y  ) ;  we  consider  approximations  for  large  scale  yielding  later.  Also,  we 
have  approximated  strains  on  the  notch  surface  by  requiring  compatibility  with 
the  homogeneous  deformation  of  the  void  interior,  eq.  (14b).  Had  we  assumed 
instead  that  the  energy  density  distribution  remains  the  same  as  for  the  linear 
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elastic  elliptical  void  problem,  eq.  (15b),  larger  values  for  emax  would  be 
obtained  in  all  cases  except  the  linear  elastic.  The  largest  difference  occurs 
in  the  elastic-perfectly  plastic  case,  for  which  eq.  (15b)  leads  to 


emax 


15  J 
8oyr 


0.50  Ey  +  2.95  Ey 


(21) 


for  loads  greater  than  the  initial  yield  value.  The  equation  is  shown  by  the 
upper  dashed  curve  in  fig.  6.  The  difference  with  eq.  (18c),  based  on  homo¬ 
geneous  deformation  of  the  void  interior,  approaches  25%  when  the  concentrated 
strain  is  large  compared  to  the  yield  strain. 

A  lower  bound  on  the  maximum  concentrated  strain  is  readily  established, 

for 


J  =  ^  Wdy  t  WUmal[>  ^  dy  =  2rt  WU^)  , 


(22a) 


with  the  latter  form  for  a  semi-circular  tip. 


For  the  elastic-perfectly  plastic  case,  the  resulting  inequality  is 


:max  *  2  (eY  +  a  r  * 


0.50  Ey  +  1.57  Ey  <5=>  - 


(22b) 


whenever  the  load  is  large  enough  for  the  right  side  of  the  inequality  to 
exceed  ey  .  This  lower  bound,  along  with  the  lower  bound  resulting  in  the 
elastic  range,  is  shown  by  the  lower  dashed  curve  in  fig.  6. 
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The  Barenblatt-Dugdale  Crack  Model 

Sharp  cracks  lead  to  strain  singularities.  The  Griffith  theory  of 
brittle  fracture  [16,17,10]  ignores  the  unrealistic  prediction  of  conditions 
at  the  crack  tip,  and  employs  an  energy  balance  to  set  the  potential  energy 
decrease  due  to  crack  extension  equal  to  the  surface  energy  of  the  newly 
created  crack  faces.  An  alternate  approach  due  to  Barenblatt  [18]  removes 
the  singularity  by  considering  a  cohesive  zone  ahead  of  the  crack,  postulating 
that  the  influence  of  molecular  or  atomic  attractions  is  representable  as  a 
restraining  stress  acting  on  the  separating  surfaces,  fig.  7.  The  restraining 
stress  o  =  o(6)  is  a  function  of  separation  distance  6  as  in  fig.  7b. 

Dugdale  [19]  employed  the  same  model  to  discuss  plane  stress  yielding  near 
cracks  in  thin  sheets.  As  further  studies  have  verified  [20],  yielding  is 
limited  to  a  narrow  slit-like  zone  ahead  of  the  crack  of  height  approximately 
equal  to  sheet  thickness  when  the  plastic  cohesive  zone  is  large  compared  to 
thickness.  Except  for  a  perturbation  near  the  tip,  yielding  then  consists  of 
slip  on  45°  planes  through  the  thickness  so  that  the  plastic  strain  is  * 
essentially  the  separation  distance  divided  by  the  sheet  thickness.  Restraining 
stresses  typical  of  plane  stress  yielding  are  shown  in  fig.  7c. 

We  may  evaluate  our  integral  J  of  eq.  (2)  by  employing  path  independence 
to  shrink  the  contour  r  down  to  the  lower  and  upper  surfaces  of  the  cohesive 
zone  as  in  fig.  7a.  Then  since  dy  =  0  on  T  , 
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[  o(<5)  d6}  dx 

•'0 


o( 6)  dfi 


(23) 


where  6^  is  the  separation  distance  at  the  crack  tip.  Thus,  if  the  crack 
configuration  is  one  of  many  for  which  J  is  known,  we  are  able  to  solve  for 
the  crack  opening  displacement  directly  from  the  force-displacement  curve. 

Now  let  us  compare  the  Griffith  theory  of  elastic  brittle  fracture  with 

the  fracture  prediction  from  the  Barenblatt  type  cohesive  force  model.  Letting 
1 fe 

6  be  the  separation  distance  in  fig.  7b  when  the  atoms  at  the  crack  tip  can 
be  considered  pulled  out  of  range  of  their  neighbors,  the  value  of  J  which 
will  just  cause  crack  extension  (or,  if  crack  extension  is  considered  as 
reversible,  will  maintain  equilibrium  at  the  current  crack  length)  is  then 

A 

fS 

J  =  a(6)  dfi  (cohesive  theory)  (24a) 

Jo 


On  the  other  hand,  the  Griffith  theory  regards  the  total  potential  energy  of 
a  cracked  body  as  P  +  2  S  &  ,  where  l  is  crack  length,  S  is  surface  energy, 
and  P  is  the  potential  energy  defined  by  the  continuum  mechanics  solution 
without  regard  to  cohesive  forces.  Determining  equilibrium  by  setting  the 
variation  in  total  potential  to  zero. 
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2  S 


(Griffith  Theory) 


(24b) 


We  conclude  that  the  Griffith  theory  is  identical  to  a  theory  of  fracture 
based  on  atomic  cohesive  forces,  regardless  of  the  force-attraction  law,  so 
long  as  the  usual  condition  is  fulfilled  that  the  cohesive  zone  be  negligible 
in  size  compared  to  characteristic  dimensions  (small  scale  yielding).  This  is 
because  the  area  under  the  force-separation  curve  is  by  definition  twice  the 
surface  energy  and,  as  we  have  noted,  J  is  equal  to  the  potential  energy 
decrease  rate  for  small  scale  yielding. 

Strain  hardening  behavior  is  also  readily  included  in  the  Dugdale  model; 
analyses  to  date  have  been  limited  to  perfect  plasticity.  For  example,  with 
linear  work  hardening 

o( 6)  =  ay  +  Ew  |  (25a) 


where  h  is  sheet  thickness,  E  a  tangent  modulus,  and  the  strain  approxi- 

w 

mated  by  6/h  ,  eq.  (23)  gives  the  crack  opening  displacement  as 


ah  2  E  J  1/2 

5t  ■  4-  {a  *  ~r u 

w  h  oy 


oy  h 
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<  1  +  2*  r- 


i) 


(25b) 


In  the  latter  form  the  plane  stress  value  of  J  for  the  crack  of  length  2a 
appropriate  for  small  scale  yielding,  has  been  used.  With  negligible  strain 
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hardening  (Ew  =  0)  this  becomes 


(25c) 


as  is  obvious  from  eq.  (23).  The  simple  picture  of  plane  stress  yielding  in 
the  Dugdale  model  is  correct  only  at  distances  from  the  crack  tip  greater  than 
the  sheet  thickness,  approximately,  with  a  truly  three  dimensional  deformation 
state  prevailing  very  near  the  crack  tip  [20].  However,  in  thin  sheets  with  the 
plastic  zone  extending  over  a  distance  several  times  the  thickness,  estimates  of 
crack  opening  displacements  may  be  expected  to  be  reasonably  accurate  as  the 
tip  displacement  is  an  integrated  effect  of  plastic  straining  along  the  entire 
length  of  the  plastic  zone.  Experimental  confirmations  of  perfectly  plastic 
predictions  for  plastic  zone  sizes  [19]  and  crack  opening  displacements  [20] 
have  been  satisfactory,  and  opening  displacements  have  proven  a  useful,  if 
imperfect,  measure  of  the  severity  of  local  deformations  in  the  formulation  of 
fracture  criteria  [21]. 

Complete  solutions  of  the  Barenblatt-Dugdale  model  have  been  given  for 
the  case  of  a  constant  cohesive  stress  aQ  acting  whenever  the  separation 
distance  6  is  non-zero.  The  resulting  opening  displacement  for  the  crack 
of  length  2a  in  a  uniform  stress  field,  fig.  4,  is  [4] 

8  (1  -  v2)  a  [  ir  am 

*t  s  — n -  log  sec  (26a) 

2 

for  plane  strain,  with  the  1  -  v  factor  removed  for  plane  stress.  But 
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since  6^  =  J/Oq  from  eq.  (23),  we  can  solve  for  J  to  find 


8  (1  -  v2)  o  2 
J  =  - -e -  log 


ir  a 

sec  (^— -) 

2  ”o 


Expanding  the  logarithm  secant  term  in  a  series  and  retaining  only  the  first 

non- vanishing  term,  as  appropriate  when  a„/°Q  is  small,  we  obtain  the  linear 

elastic  expression  for  J  of  eq.  (9)  and  verify  our  procedure  of  using  the 

linear  elastic  expressions  in  cases  of  small  scale  yielding.  Solutions  of  the 

Barenblatt-Dugdale  model  with  constant  cohesive  stress  have  been  given  for 

other  configurations  of  single  cracks  in  infinite  bodies  [4],  including  wedge 

force  loadings,  and  for  an  infinite  periodic  array  of  collinear  cracks  [22], 

an  approximation  for  single  edge,  double  edge,  and  central  cracks  in  finite 

width  plates.  Crack  opening  displacements  from  these  solutions  may  be  used 

to  compute  J  ,  as  in  eqs.  (26).  In  view  of  the  physical  interpretation 

given  J  as  an  energy  variation  rate,  we  do  not  expect  resulting  expressions 

to  be  too  tied  to  the  model.  Thus,  after  a  judicious  choice  of  to  account 

for  possible  hydrostatic  stress  elevations  of  plastic  limit  loads  as  in  the 

double  edge  notched  plate  [23],  Barenblatt-Dugdale  estimates  of  J  may  be 

used  as  an  approximation  in  cases  of  large  scale  yielding  for  computations 

such  as  those  of  the  last  section.  Fig.  8  shows  a  comparison  of  the  linear 

elastic  (eq.  9)  and  constant  cohesive  stress  Barenblatt-Dugdale  (eq.  26b) 

expressions  for  J  for  the  configuration  of  fig.  4.  The  plot  is  in  dimension- 

2  2 

less  form,  with  E  J/[ir(l  -  v  )  oQ  ]  shown  as  a  function  of  the  applied  stress 
ratio  a  /o.  .  The  linear  elastic  (small  scale  yielding)  result  is  a  parabola 

00  (J 
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with  the  dimensionless  J  equal  to  the  square  of  the  stress  ratio.  The 
closeness  of  the  two  curves  at  low  stress  levels  indicates  the  great  use¬ 
fulness  of  the  small  scale  yielding  approximation. 


Because  it  will  be  useful  for  a  later  comparison  with  plane  strain 
plasticity,  we  record  here  the  expression  for  the  length  of  the  cohesive  zone 
in  the  case  of  constant  cohesive  stress  and  the  crack  of  length  2a  in 

a  uniform  stress  field.  Letting  R  be  the  cohesive  zone  length,  the  result 
is  [4] 


R 


it  a 

sec  (y-  °°) 
2  °0 


(26c) 


The  latter  form  is  the  first  term  in  a  Taylor  expansion.  Again  we  verify 
that  the  elastic  stress  intensity  factor  (eq.  9)  controls  inelastic  behavior 
in  the  small  scale  yielding  range.  Setting  =  Oy  (the  yield  stress  in 
tension),  R  is  the  plastic  zone  size  in  plane  stress  predicted  by  the  Dug- 
dale  model.  The  model  is  valid  when  the  computed  value  of  R  is  several 
times  the  sheet  thickness. 

Other  known  elastic-plastic  solutions  for  cracks  and  notches  may  be  used 
to  estimate  J  for  cases  of  large  scale  yielding.  In  particular  the  many 
available  perfectly  plastic  [5,7,8]  and  strain  hardening  [2,6,8]  solutions 
for  cracks  in  anti-plane  stress  fields  may  be  put  to  use.  Resulting  detailed 
evaluations  of  J  are  deferred  in  the  interest  of  brevity.  It  is  simply  noted 
here  that  in  all  cases  the  validity  of  our  small  scale  yielding  result  (last 
term  of  eq.  10)  is  confirmed,  and  that  in  the  perfectly  plastic  case  the  large 
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scale  yielding  range  prediction  of  the  deviation  of  J  from  the  small  scale 
yielding  formula  is  not  significantly  different  from  that  of  the  Barenblatt- 
Dugdale  model. 

Plane  Strain  Yielding  Near  a  Crack 

An  approximate  treatment  of  perfectly  plastic  plane  strain  yielding  near 
a  crack  may  be  based  on  the  slip  line  theory  [15]  and  our  considerations  of 
path  independence.  Slip  line  theory  is  really  a  valid  consequence  of  a  Tresca 
or  Mises  yield  criterion  only  when  elastic  straining  is  incompressible  or  when 
plastic  strains  greatly  exceed  initial  yield  strains.  However,  we  shall  see 
that  it  is  appropriate  for  our  present  purposes.  Suppose  tentatively  that  the 
plastic  region  completely  surrounds  the  crack  tip  as  in  fig.  9.  The  exact 
shape  of  the  elastic-plastic  boundary  is  unknown,  as  indicated  by  the  question 
mark.  The  traction  free  boundary  condition  on  the  crack  surface  determines  a 
constant  stress  state  in  the  largest  isosceles  right  triangle  (A  in  fig.  9) 
which  may  be  fit  into  the  plastic  zone,  and  in  this  region 

o  =2tv  ,  o  =o  =0  (region  A)  (27a) 

xx  Y  *  yy  xy 

where  Ty  is  the  yield  stress  in  pure  shear.  We  choose  tension  instead  of 
compression  for  a  since  the  latter  choice  would  result  in  biaxial  com¬ 
pression  ahead  of  the  crack.  Any  slip  line  emanating  from  the  crack  surface 
and  finding  its  way  to  the  x  axis  in  front  of  the  crack  must  swing  through 
an  angle  of  ir/2  .  Thus  a  hydrostatic  stress  elevation  of  2  Xy(ir/2)  must 
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result  ahead  of  the  crack  [15].  Constant  stresses  on  this  line  determines 
a  constant  stress  state  in  a  diamond  shaped  region  (B  in  fig.  9),  and  there 

a  =  ttt„  ,  o  =  (2  +  ir)  tv  ,  a  =0  (region  B)  .  (27b) 

xx  Y  *  yy  Y  ’  xy  & 

A  centered  fan  (C  in  fig.  9)  must  join  two  such  regions  of  constant  stress 
[15,24]  and  stresses  in  the  fan  are 

°rr  *  o60  =  ^  +  Ty  ~  2ty6  »  are  =  ty  (region  C).  (27c) 

The  resulting  slip  line  field  is  familiar  in  the  limit  analysis  of  rigid 
indentors  [15,24]  and  double  edge  notched  thick  plates  [23,25]. 

Large  strains  can  occur  only  when  slip  lines  focus,  as  in  the  centered 
fan,  but  not  in  constant  state  regions  (A  and  B)  unless  strains  are  uniformly 
large  along  the  straight  slip  lines  and  thus  out  toward  the  elastic-plastic 
boundary  (where  strains  must  be  small).  One  then  anticipates  strains  on  the 
order  of  initial  yield  values  in  the  constant  state  regions.  As  we  shall  see, 
only  the  strain  singularity  in  the  centered  fan  enters  the  results  we  are 
heading  toward.  Assuming  elastic  incompressibility  and  agreement  between 
principal  shear  stress  and  strain  directions,  e.^  =  =0  in  the  centered 

fan.  This  implies  that  displacement  components  are  there  representable  in  the 
form  [15,24] 

u,  =  f'(6)  ,  u.  =  -f( 6)  +  g(r)  (28a) 

V  V 
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Note  that  in  a  proper  incremental  theory,  these  equations  apply  to  velocities 
instead  of  displacements.  The  shear  strain  in  the  centered  fan  is 

given  by 


e 


r6 


1  (I 

2  vr 


3ue 

3r 


f"(6)  +  f(0)  -  g(r)  .  1  _',_s 

- - -  +  2g(r)  * 

(28b) 


Now  consider  the  path  independent  integral  J  ,  eq.  (2).  Taking  T 
to  be  a  circle  of  radius  r  centered  at  the  crack  tip  and  expressing  all 
quantities  in  polar  component  form,  there  results 

f+ir 

J  =  r  {W  cos  6  -  c  [e  cos  6  -  (e  „  -  <o)  sin  0] 

l  rr  rr  ru 

*  -TT 

-  a  .  [(e  Q  +  w)  cos  0  -  eQO  sin  0]}  d0  .  (29a) 

Tv  To  oo 


Here  u  is  the  rotation,  measured  positive  contraclockwise : 


b) 


1  f3ue 

2  ^3r 


u. 


i  3u 

1  _ r 

r  36 


(29b) 


We  evaluate  the  integral  by  letting  r  -*•  0  (as  permitted  by  path  independ¬ 
ence).  Since  strains  are  non-singular  in  the  constant  state  regions,  only 
the  portion  of  the  integral  taken  over  the  centered  fans  (ir/4  <  |  0 1  <  3ir/4) 
contributes  in  the  limit.  We  now  express  the  energy  density,  strains,  and 
rotation  in  a  form  appropriate  when  r  ■*>  0  .  First  note  that  arbitrarily  we 
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may  set  f(ir/4)  =  0  .  Then  g(r)  represents  the  0  direction  displacement 
of  the  straight  line  boundary  between  regions  B  and  C  in  fig.  9  .  Since 
strains  are  bounded  in  the  constant  state  region,  the  displacement  is  zero  at 
r  =  0  and  0  =  *ir/4  .  Thus 

g(0)  =  0  and  g’(0)  =  r1^iQ  ,  (30a) 

and  expressions  for  displacements,  strains,  and  rotation  very  near  the  crack 
tip  become 


u  =  f'(0)  ,  ue  =  -f  ( 0)  (fCir/4)  =  f'(ir/4)  =  0) 


_  f  (0)  +  f ( 0 ) 

err  “  E00  "  0  »  er0  “  ““  “  2  r 


(30b) 


The  energy  density  for  an  elastically  incompressible  perfectly  plastic 
material  is 


12  1 
W  =  —  Gy  (y  <  yy)  ,  W  =  T^y  -  tyYy  (Y  *  Yy) 


Here  y^  =  Ty/G  is  the  "engineering"  yield  strain  in  shear  and 


Ola) 


Y  = 


(e  -  £00)  +  4e 

rr  00  r0 


1/2 


(31b) 


is  the  principal  shear  strain.  We  may  see  that  the  expression  given  for  W 
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when  Y  >  Yy  is  correct  in  the  sense  that  any  set  of  stresses  satisfying 

dW  =  o.,de..  also  satisfies  the  Mises  or  Tresca  yield  condition  and  gives  a 
13  iD 

common  ratio  of  deviatoric  stress  components  to  deviatoric  strain  components . 
Thus  as  r  **■  0  in  the  centered  fans 


W 


f"( e)  +  f(e) 


r 


(31c) 


Letting  r  -*■  0 


in  eq.  (29a),  the  result  for  J  then  becomes 


J  =  2t, 


3ir/4 

tt/4 


[f  (0)  +  f ( 0)3  [cos  0  +  (1  +  -^ - 20)  sin  0]  d0 


(32) 


Regarding  J  as  known,  we  get  an  approximate  solution  by  picking  a  reasonable 
form  for  f(0)  containing  an  unknown  constant  and  determine  the  constant  by 
evaluating  the  integral.  Before  a  bit  of  further  analysis  to  identify  reason¬ 
able  choices,  we  review  our  assumptions.  Elastic  compressibility,  rather  than 
our  assumed  incompressibility,  appears  unimportant,  at  least  in  eq.  (32),  as 
only  the  strain  singularity  contributes  and  slip  line  theory  is  valid  when 
strains  are  large  compared  to  initial  yield  strains.  The  more  serious  assumption 
is  that  the  stress  field  given  by  eqs.  (27)  is  assumed  to  actually  occur  near  the 
crack  tip.  As  indicated,  this  stress  field  is  correct  if  the  plastic  zone  com¬ 
pletely  surrounds  the  crack  tip.  Some  as  yet  unpublished  etching  experiments  by 
Prof.  F.  A.  McClintock,  discussed  in  [8],  reveal  a  plastic  zone  that  does  sur¬ 
round  the  tip.  But  other  etching  experiments  [26]  are  inconclusive  on  whether 
yielding  actually  occurs  over  any  sizable  region  directly  ahead  and  behind  the 
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crack  tip.  It  is  worth  noting,  however,  that  the  elastic-plastic  boundary 
could  cut  sharply  into  the  crack  tip  in  regions  A  and  B  in  fig.  9 
without  effecting  the  validity  of  our  stress  field  very  near  the  tip.  In 
particular,  the  stresses  given  by  eqs.  (27a  and  b)  could  be  valid  right  at 
the  crack  tip  in  elastic  material,  with  a  steady  decay  of  the  resolved  shear 
stress  away  from  the  yield  value  at  non-zero  distances  from  the  tip  in  regions 
A  and  B  .  The  centered  fan  could  then  still  result,  and  eq.  (32)  remain 
valid.  It  is  especially  interesting  that  a  result  of  just  this  sort  occurs 
for  the  perfectly  plastic  anti-plane  strain  crack  problem.  In  that  case  a 
construction  of  the  stress  field  on  the  assumption  of  a  plastic  zone  sur¬ 
rounding  the  tip  leads  to  a  fan  of  anti-plane  shear  lines  centered  at  the 
crack  tip  in  the  region  x  >  0  of  fig.  9  (oa  =  rv)  ,  and  constant  stress 
states  adjoining  the  crack  surfaces  for  x  <  0  (o  =  *!„)  .  But  exact 

XZ  Y 

solutions  [2,5,7]  reveal  an  elastic-plastic  boundary  extending  ahead  of  the 
crack  and  cutting  into  the  crack  tip  tangentially  to  the  boundary  between  the 
centered  fan  and  constant  stress  regions.  The  true  plastic  zone  encloses  only 
points  in  the  centered  fan  region;  in  what  was  apparently  a  constant  stress 
plastic  region,  there  is  elastic  behavior  with  a  steady  decay  of  the  resolved 
shear  stress  away  from  the  yield  value  as  one  moves  away  from  the  crack  tip. 

Let  us  define  a  function  R(0)  by 

Yy  R(6)  =  f"(e)  +  f ( 0 )  .  (33a) 

Then  from  eqs.  (30b, 31b),  the  "engineering"  shear  strain  near  the  crack  tip 
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in  the  fan  is  given  by 


Y  =  Yy 


R(e) 

r 


(33b) 


Note  from  eq.  (28b)  that  this  equation  would  apply  throughout  the  centered 
fan  if  g(r)  were  a  linear  function  of  r  (as  would  be  the  case  if  the  con¬ 
stant  state  region  B  in  fig.  9  underwent  a  homogeneous  strain).  Thus  R(0) 
is  an  approximate  indication  of  the  distance  to  the  elastic-plastic  boundary, 
and  we  can  estimate  the  plastic  zone  size  by  writing  eq.  (32)  in  the  form 

f3ir/4 

J  =  2t  y  R(6)  [cos  6  +  (1  +  ^  -  26)  sin  6]  d6  .  (33c) 

•'ir/4 

For  a  given  choice  of  R(0)  ,  displacements  near  the  crack  tip  are  computed 
from  eqs.  (30b)  by  solving  the  differential  equation  (33a)  under  initial  con- 

f 

ditions  f(ir/4)  =  f  ( tt/4 )  =  0  .  It  is  also  useful  to  write  eq.  (32)  in  terms 
of  cartesian  components  of  displacements : 

u  =  u  sin  0  +  un  cos  0 

y  r  0 

u  =  u  cos  0  -  u„  sin  6  (34a) 

x  r  0 


Differentiating  with  respect  to  0  and  using  eqs.  (30b),  in  the  fan 


du  „ 

=  [f  (6)  +  f(0)]  sin  0  =  yy  R(6)  sin  0 


du 


x 


d0 


=  [f  (0)  +  f(6)]  cos  0  = 


Yy  R(0)  cos  0 


(34b) 
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The  ordinary  differentiation  symbol  is  used  as  displacements  depend  only 
on  0  very  near  the  tip.  We  can  then  estimate  crack  opening  displacements 
by  writing  eq.  (32)  in  the  form 

/■3ir/4  du 

J  =  2tY  Jir/4  *-sin  6  +  (1  +  2~  ~  2e)]  d6  *  (34c) 

Let  us  first  consider  the  crack  opening  displacement  6^  =  2u^  at  0  =  3ir/4 
Assume  as  an  approximation  that  du^/d0  is  symmetric  about  0  =  ir/2  .  From 
eq.  (34b),  this  is  equivalent  to  assuming  symmetry  of  R(0)  about  the  vertical: 
R(ir/2  +  0)  =  R(ir/2  -  0).  The  bracketed  term  in  eq.  (34c)  may  be  split  into  a 
symmetric  part  equal  to  (1  +  ir/2)  and  an  anti-symmetric  part.  The  anti¬ 
symmetric  contributes  zero  to  the  integral  and  thus 

1 3ir/4  du 

J  <=  2tv  (1  +  tt/2)  d0  =  (1  +  ir/2)  tv6 

Y  Jir/4  d0  Y  * 

2tt(1  -  v2)o2a 

x  _  Z  U  _  00  /  nr  \ 

t  ~  2  +  TT  Ty  =  (2  +  Tl)  E  Ty 

Here,  as  throughout  the  paper,  the  latter  form  applies  for  small  scale  yielding 
near  a  crack  of  length  2a  in  a  uniform  stress  field  om  .  Comparing  with  the 
Dugdale  plane  stress  value  (eq.  25c)  for  the  same  remotely  applied  stress  and 
crack  length,  and  for  \>  =  0.3  ,  the  plane  strain  crack  opening  displacement  is 
61%  of  the  plane  stress  value  for  a  Mises  yield  condition  (cry  =  /3  Ty)  and 
70%  for  a  Tresca  yield  condition  (oy  =  2ty)  .  Note  that  if  R(0)  is  symmetric 
about  0  =  ir/2  ,  as  we  have  just  assumed,  then  an  integration  of  the  latter  of 
eqs.  (34b)  gives  u^  =  0  at  0  =  3ir/4  .  This  means  that  the  crack  opens  up 
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into  a  rectangular  shaped  tip,  much  as  observed  by  Laird  and  Smith  [26]. 

An  integration  by  parts  in  eq.  (34c)  expresses  J  in  terms  of  an  integral 
of  Uy  times  a  function  of  6  .  From  eq.  (34b)  du^/de  is  non-negative 
since  R(0)  is  non-negative,  so  that  u^  $  6^/2  .  Applying  this  inequality 
to  the  expression  for  J  (after  integration  by  parts),  one  obtains  a  lower 
bound  on  the  crack  opening  displacement: 


6 


t 


% 


1  J__ 

2  +  TT  Ty 


(36) 


Our  reason  for  concern  with  crack  opening  displacements  is  shown  in 
fig.  10.  The  slip  line  field  of  fig.  9  suggests  no  intense  strain  concentration 
ahead  of  the  crack.  But  when  the  slip  line  field  is  drawn  for  a  crack  blunted 
by  plastic  deformation,  a  very  different  picture  results  on  a  small  scale  of 
the  order  of  6  .  The  fan  C  becomes  non-centered  and  its  straight  slip  lines 

focus  into  a  region  D  of  intense  deformation.  For  simplicity  of  illustration, 
the  crack  tip  has  been  drawn  as  a  semi-circle  in  fig.  10  and  the  associated 
exponential  spiral  slip  line  field  [15,24]  extends  a  distance  of  1.9  6^. 
ahead  of  the  blunted  tip.  From  estimates  of  the  plastic  zone  size  to  be  given 
shortly,  we  shall  see  that  6^  is  of  the  order  of  the  initial  yield  strain 
times  a  linear  dimension  of  the  plastic  zone,  so  that  the  intense  deformation 
region  is  extremely  small  and  fig.  10  is  essentially  fig.  9  with  a  magnification 
in  linear  dimensions  of  the  order  of  one  over  the  initial  yield  strain.  Since 
the  blunted  region  is  small,  an  effective  procedure  would  be  to  perform  an 
incremental  analysis  of  blunting  by  regarding  the  constant  displacement  rate 
along  each  straight  slip  line  of  the  noncentered  fan  to  be  given  by  the  rate 
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of  increase  of  u  =  u  (0)  of  our  present  analysis  (where  8  is  now  the 
r  r 

inclination  of  a  given  straight  slip  line  and  is  identical  to  the  polar 
coordinate  0  at  distances  from  the  tip  large  compared  to  6^)  . 

Calling  R(0)  the  distance  from  the  crack  tip  to  the  elastic-plastic 
boundary,  as  an  approximation,  let  us  assume  that 

R(0)  =  R  (a  constant)  (37a) 


so  that  the  boundary  is  an  arc  of  circle  of  radius  R  in  the  centered  fan. 
Substitution  into  eq.  (33c)  estimates  the  plastic  zone  dimension  as 


R 


GJ _ 

/2  (2  +  ir) 


it(1  -  v) 

2  S2  (2  +  ir) 


(37b) 


We  shall  see  that  this  is  actually  the  smallest  possible  value  which  the 
maximum  value  of  R( 0 )  may  have .  Near  crack  tip  displacements  in  the  fan 
associated  with  this  choice  of  R(0)  are 

ur  =  YyR  sin  (0  -  tt/4) 

u  =  -YyR  Cl  -  cos  (0  -  tt/4)]  (37c) 

Another  choice  for  R(0)  , 

R(  0)  =  R  sin  (28  -  ir/2)  (R  a  constant)  (38a) 
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gives  an  elastic-plastic  boundary  which  cuts  into  the  crack  tip  along  radial 
boundaries  of  the  centered  fan  (R(tt/4)  =  R(3ir/4)  =  0)  as  discussed  earlier. 
In  this  case  eq.  (33c)  estimates  the  maximum  plastic  zone  dimension  as 


2/2  (2  +  it)  t. 


3tt(  1  -  v) 
4/2  (2  +  it) 


(38b) 


50%  greater  than  the  estimate  of  eq.  (37b)  above.  Associated  near  tip  dis- 
placements  are 

2yyR 

u  =  — 5—  [cos  (0  -  tt/4)  -  cos  (20  -  it/2)] 

r  o 


yyR 

u.  =  -  —5—  [2  sin  (0  -  tt/4 )  -  sin  (20  -  tt/2)]  .  (38c) 

U  o 


Comparing  with  the  displacements  of  eq.  (37c),  both  give  the  crack  opening 
displacement  of  eq.  (35),  as  they  must  since  R(0)  is  chosen  symmetrical 
about  0  =  tt/2  in  both  cases,  u^  is  the  input  to  the  large  deformation 
region  in  an  analysis  of  blunting,  and  the  two  predictions  of  up  agree 
identically  for  0  =  tt/4  ,  tt/2  ,  and  3ir/4  .  The  greatest  difference  occurring 
for  intermediate  angles  is  about  18%  of  the  maximum  value  of  u^  =  6^/(2/2) 
occurring  at  0  =  3tt/4  ,  so  predictions  of  blunting  appear  insensitive  to  the 
choice  of  R(0)  . 

Comparing  the  two  estimates  of  plastic  zone  size  (eqs.  37b  and  38b)  with 
the  plane  stress  zone  size  (eq.  26c)  predicted  by  the  Dugdale  model,  for 
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v  =  0.3  eq.  (37b)  predicts  a  plane  strain  R  which  is  37%  of  the  plane 
stress  R  for  a  Mises  material  and  49%  for  a  Tresca  material.  Corres¬ 
ponding  figures  for  eq.  (38b)  are  55%  and  73%  .  Observations  [27]  by 
etching  suggest  a  figure  in  the  neighborhood  of  50%  .  Eq.  (33c)  gives  a 
lower  bound  on  the  maximum  value  of  R(0)  occurring  in  the  fan.  Inserting 
the  inequality  R(0)  S  ^max  »  one  obtains 


R 


max 


GJ 

/2  (2  +  ir) 


ir(l  -  v) 
2/2  (2  +  it) 


(39) 


The  lower  bound  is  identical  to  the  value  of  R  given  by  eq.  (37b)  . 
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NOTCH  / 


fig.  1.  Flat  surfaced  notch  in  two-dimensional 
deformation  field  (all  stresses  depend  only  on 
x  and  y).  r  is  any  curve  surrounding  the  notch 
tip;  r  denotes  the  curved  notch  tip. 
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fig.  2.  Two  special  configurations  for  which  the  path  independent  integral 
J  is  readily  evaluated  on  the  dashed  line  paths  T  shown.  Infinite  strips 
with  semi-infinite  notches,  (a)  constant  displacements  imposed  by  clamping 
boundaries,  and  (b)  pure  bending  of  beam-like  arms. 


fig.  3.  (a)  Small  scale  yielding  near  a  narrow  notch  or  crack  in  an  elastic- 

plastic  material,  (b)  The  actual  configuration  is  replaced  by  a  semi-infinite 
notch  or  crack  in  an  infinite  body,  actual  boundary  conditions  are  replaced  by 
the  requirement  of  an  asymptotic  approach  to  the  linear  elastic  crack  tip 
singularity  stress  field. 


fig.  4.  Narrow  notch  or  crack  of 
length  2a  in  infinite  body. 
Uniform  remote  stress  o 


fig.  5.  Coordinates  employed  in 
description  of  notch  surface. 

4>  is  tangent  angle  and  rt(<{>) 
is  radius  of  curvature. 
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ELASTIC  -  PERFECTLY  PLASTIC  STRAIN  CONCENTRATION 
FLAT  SURFACED  NOTCH,  SEMI-CIRCULAR  TIP. 


fig.  6.  Strain  concentration  by  flat  surfaced  notch  with  semi-circular  tip  of 
radius  r  in  perfectly  plastic  material. 


(c)  PLANE  STRESS  YIELDING 


fig.  7.  Dugdale-Barenblatt  crack  model,  (a)  Cohesive  zone  at  crack  tip  with 
restraining  stress  dependent  on  separation  distance,  (b)  Force-displacement 
relation  for  atomic  attraction  in  elastic  brittle  fracture,  and  (c)  for  plane 
stress  plastic  yielding  in  thin  sheet. 


PATH  INDEPENDENT  INTEGRAL,  J  (DIMENSIONLESS) 


EVALUATION  OF  PATH  INDEPENDENT  INTEGRAL  J 
FOR  CRACK  OF  LENGTH  2a  IN  INFINITE  BODY  UNDER  STRESS  CT^. 


fig.  8.  Comparison  of  small  scale  yielding  result  for  J  with  exact  value 
from  Barenblatt  -  Dugdale  crack  model  (can  be  used  to  approximate  J  in  the 
large  scale  yielding  range  for  other  problems  of  cracks  or  narrow  notches). 
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fig.  9.  Perfectly  plastic  plane  strain  slip  line  field  at  a  crack  tip. 
Constant  stress  regions  A  and  B  joined  by  centered  fan  C  .  Position 
of  the  elastic-plastic  boundary  is  unknown;  it  may  cut  in  toward  the  crack 
tip  in  regions  A  and  B  (see  discussion). 


fig.  10.  Crack  tip  blunting  creates  a  small  region  D  of  intense  deforma¬ 
tion  ahead  of  the  crack.  This  is  fig.  9  magnified  in  linear  dimensions 
by  a  large  factor  of  the  order  of  one  over  the  initial  yield  strain. 


